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Abstract 

Using the two-loop approximation of chiral perturbation theory, we calculate the mo- 
mentum and density dependent isovector nuclear spin-orbit strength Vu(p, kf). This quan- 
tity is derived from the spin-dependent part of the interaction energy T, sp i n = ■ (q x 
p)[Ui s (p, kf) — Vi s (p, kf)r^5\ of a nucleon scattering off weakly inhomogeneous isospin- 
asymmetric nuclear matter. We find that iterated l7r-exchange generates at saturation 
density, kfo = 272.7 MeV, an isovector nuclear spin-orbit strength at p = of Vi s (0, kfo) ~ 
50MeVfm 2 . This value is about 1.4 times the analogous isoscalar nuclear spin-orbit 
strength J7/ s (0, /c/o) ~ 35MeVfm 2 generated by the same two-pion exchange diagrams. 
We also calculate several relativistic 1/M-corrections to the isoscalar nuclear spin-orbit 
strength. In particular, we evaluate the contributions from irreducible two-pion exchange 
to Ui s (p,kf). The effects of the three-body diagrams constructed from the Weinberg- 
Tomozawa 7T7riViV-contact vertex on the isoscalar nuclear spin-orbit strength are com- 
puted. We find that such relativistic 1/M-corrections are less than 20% of the isoscalar 
nuclear spin-orbit strength generated by iterated one-pion-exchange, in accordance with 
the expectation from chiral power counting. 

PACS: 12.38.Bx, 21.65.+f, 24.10.Cn 
Keywords: Effective field theory at finite density, Isoscalar and isovector nuclear spin-orbit 
interaction. 



1 Introduction 



The introduction of the spin-orbit term into the nuclear single-particle Hamiltonian has been 
most decisive for the success of the nuclear shell model [1]. Only with a very strong and 
attractive spin-orbit potential one is, for example, able to explain the observed sequence of so- 
called magic numbers {2, 8, 20, 28, 50, 82, 126, . . .}. The dynamical origin of the strong nuclear 
spin-orbit force has not been fully resolved even up to date. The analogy with the spin-orbit 
interaction in atomic physics gave the hint that it could be a relativistic effect. This idea has 
lead to the construction of the relativistic (scalar- vector) mean-field model [2]. In this model 
the spin-independent nuclear potential (of approximate depth —50 MeV) results from an almost 
complete cancelation of a very strong attraction generated by scalar (a-meson) exchange and 
a nearly equally strong repulsion generated by vector (tu-meson) exchange. The corresponding 
spin-orbit term (obtained by a non-relativistic reduction of the nucleon's Dirac-Hamiltonian) 
comes out proportional to the coherent sum of the very large scalar and vector mean-fields. 
In this sense, the relativistic mean-field model gives a simple and natural explanation of the 
basic features of the nuclear shell model potential. Refinements of relativistic mean-field models 



which include additional non-linear couplings of the scalar and vector fields or explicitly density- 
dependent point couplings of nucleons are nowadays widely and successfully used for nuclear 
structure calculations [3, 4, 5]. 

The nuclear spin-orbit potential arises generally as a many-body effect from the underlying 
spin-orbit term in the (free) nucleon-nucleon scattering amplitude. The calculation of the tree 
level diagrams with one scalar-meson or one vector-meson exchange between nucleons gives 
indeed a spin-orbit term in the NN T-matrix proportional to 1/M 2 , with M denoting the nucleon 
mass. The nuclear spin-orbit potential corresponding to scalar and vector meson exchange 
is therefore obviously a truly relativistic effect. However, the quadratic reciprocal scaling of 
the spin-orbit NN-amplitude with the nucleon mass M is not universal, and it changes if one 
considers the exchange of two mesons between nucleons, i.e. loop diagrams. For example, 
irreducible two-pion exchange gives rise to a spin-orbit term in the NN T-matrix proportional 
to 1/M (see eqs. (22,23) in ref.[6]) and iterated one-pion exchange produces a spin-orbit term in 
the NN T-matrix which even scales linearly with the nucleon mass M (see eq.(33) in ref.[6]). 

In a recent work [7] we have used the systematic framework of chiral perturbation theory 
to calculate the nuclear spin-orbit interaction generated by one- and two-pion exchange. The 
momentum and density dependent nuclear spin-orbit strength Ui s (p, kf) is derived from the spin- 
dependent part of the interaction energy £ Sf ,j„ = | a • (q x p) Ui s (p, kf) of a nucleon scattering 
off weakly inhomogeneous isospin-symmetric nuclear matter. When working to linear order 
in small spatial gradients (which corresponds to linear order in the momentum transfer q) no 
detailed knowledge about the density distribution of nuclear matter is required. In fact, after 
isolating the important proportionality factor q from the expressions of self-energy-type diagrams 
the nuclear spin-orbit strength Ui s (p, kf) can be finally computed in the limit of homogeneous 
isospin-symmetric nuclear matter (characterized by its Fermi momentum kf). It was found 
in ref.[7] that the Hartree and Fock diagrams of iterated one-pion exchange generate a spin- 
orbit strength at p = of Ui s (0,kfo) = 34.7MeVfm 2 which is in perfect agreement with the 
empirical value Ui s ~ 35MeVfm 2 [9, 10] used in shell-model calculation of nuclei. This novel 
spin-orbit strength is neither of relativistic nor of short range origin. It is in fact linearly 
proportional to the nucleon mass M and its inherent range is the pion Compton wavelength 
m^ 1 = 1.46 fm. A strong p-dependence of the nuclear spin-orbit strength Ui s (p, kfo) arising from 
iterated l7r-exchange has however been observed in ref.[7]. A further property of the spin-orbit 
Hamiltonian emerging from that diagrammatic calculation was that it has (in coordinate space) 
terms proportional to V/(r) as well as terms proportional to f(r) V/(r) (with f(r) = p(r)/p(0) 
the normalized radial density profile) which get weighted differently at the surface of a finite 
nucleus. All such (unusual) features which go beyond the simple shell model parameterization 
of the nuclear spin-orbit Hamiltonian leave questions about the ultimate relevance of the spin- 
orbit interaction generated by 27r-exchange for a finite nucleus. Nuclear structure calculation 
which use the results of ref . [7] as input are necessary in order to clarify the role of the nuclear 
spin-orbit interaction generated by 27r-exchange. 

The purpose of the present work is to calculate using the same framework as in ref. [7] the 
isovector spin-orbit strength Vi s (p, kf) generated by iterated l7r-exchange in order to reveal 
the underlying isospin dependence. Furthermore, we calculate here several relativistic 1/M- 
corrections to the isoscalar nuclear spin-orbit strength Ui s (p,kf). In particular, we evaluate 
the contributions to Ui s (p,kf) arising from irreducible two-pion exchange. We compute the 
isoscalar nuclear spin-orbit strength generated by certain three-body diagrams constructed from 
the leading order isovector chiral 7T7riViV-contact vertex (so-called Weinberg- Tomozawa vertex). 
The relativistic 1/M-corrections to the dominant Hartree diagrams of iterated l7r-exchange are 
also calculated. These results shed some light on the convergence of the small momentum 
expansion underlying the whole approach. 



2 General considerations about the nuclear spin-orbit 
term 



Let us begin with recalling the spin-orbit Hamiltonian of the nuclear shell model [9, 10] which 
is generally written in the form: 



n la = (u la -v l8 T 3 s) 



G 



df(r) 



2r dr 



(1) 



Here, a is the conventional Pauli spin-vector and i = f x p is the orbital angular momentum 
of a nucleon. f(r) denotes the normalized radial density distribution of a nucleus, typically 
parameterized by a Saxon- Woods function. 5 — (N — Z) / \N + Z) is the isospin asymmetry 
parameter and r 3 — > ±1 for a proton or a neutron. The empirical value of the isoscalar nuclear 
spin-orbit strength is [/j s ~ 35MeVfm 2 [9, 10]. The ratio of the isovector and isoscalar nuclear 
spin-orbit strengths Vi s /Ui s ~ —0.6 is commonly assumed to be equal to the ratio of the spin- 
independent isovector and isoscalar average nuclear potentials [9, 10]. 

We wish to calculate the nuclear spin-orbit strengths C/j s and V\ s (or appropriate generaliza- 
tion of them) in the systematic framework of in-medium chiral perturbation theory [8]. The first 
observation one makes is that the spin-orbit interaction vanishes identically in infinite homoge- 
neous nuclear matter since there is no preferred center in this system in order to define an orbital 
angular momentum. Therefore one has to generalize the calculation of the single-particle poten- 
tial in ref.[8] to (at least) weakly inhomogeneous nuclear matter. The relevant quantity in order 
to extract the isoscalar and isovector nuclear spin-orbit strengths is the spin-dependent part 
of the interaction energy of a nucleon scattering off weakly inhomogeneous isospin-asymmetric 
nuclear matter from initial momentum p — q/2to final momentum p + q/2, which reads: 



Zspin = t -a-{qxp) {U ls {p, k f ) - V ls (p, k f ) r 3 5 + 0{5 2 )} 



(2) 



The (small) momentum transfer q is provided by the Fourier-components of the inhomoge- 
neous nuclear matter distribution. The density form factor $(<f) = / d 3 r e~' l ^ r f(r) plays the 
role of a probability distribution of these Fourier-components. With the help of this rela- 
tionship the spin-orbit Hamiltonian H\ s in eq.(l) becomes equal to the Fourier-transform of 
the product of the density form factor $(<f) and the spin-dependent interaction energy S sp j„: 
Hi, = (27r)~ 3 /rf 3 ge^$(g)S, 
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Fig.l: The double line symbolizes the medium insertion for a weakly inhomogeneous many- 
nucleon system T(pj,qj) defined by eqs.(3,4)- Right: The ln-exchange Fock graph. 



Consistent with the assumption of a weakly inhomogeneous nuclear matter distribution we 
keep in S sp j„ only linear terms in q corresponding to small spatial gradients of the density 
profile f(r). Practically, this means that after isolating the proportionality factor q in an 
explicit calculation one can take the limit of homogeneous isospin-asymmetric nuclear matter 
characterized by its proton and neutron Fermi momenta k Pj7l = kf(l =F 5) 1//3 - The expansion up 
to linear order in the isospin asymmetry parameter S gives then the momentum and density 
dependent isoscalar and isovector nuclear spin-orbit strengths Ui s (p, kf) and Vi s (p, kf). Possible 
higher powers in q contributing to S sp i„ have no direct correspondence to the (standard) spin- 
orbit Hamiltonian 7ii s in eq.(l). 

In ref. [8] the calculation of the spin-independent single-particle potential in homogeneous 
isospin-symmetric nuclear matter has been organized in the number of so-called medium in- 
sertions. The medium insertion for a (non-relativistic) many-fermion system is generally con- 
structed from the sum over the occupied energy eigenstates as [11, 12]: 

V[pj,qj) = f d 3 n f d 3 r 2 £ M^WEin)^^ 1 -^ e-i^ +r ^ . (3) 

J J E£occ 

The double line in the left picture of Fig. 1 symbolized this medium insertion together with the 
assignment of the out- and in-going nucleon momenta pj ± qj/2. The momentum transfer cfj is 
provided by the Fourier-components of the inhomogeneous matter distribution. Semi-classical 
expansions [11, 12] give for a weakly inhomogeneous and spin-saturated many-nucleon system: 

r(^,^) = |^.|) + 1^^^(A;„- (4) 

with the density form factor introduced after eq.(2). The subleading O(qj) term in eq.(4) 

will never come into play in the calculation of the spin-dependent interaction energy S sp j„ to 
linear order in q. 



3 Diagrammatic calculation of the isovector nuclear 
spin-orbit strength 
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Fig. 2: The iterated lir-exchange Hartree graphs. These diagrams are labeled (a), (b), (c), 
(d) from left to right. 

In this section we present analytical results for the isovector nuclear spin-orbit strength 
Vi s (p, kf) as given by chiral one- and two-pion exchange. We start with the l7r-exchange Fock 1 

1 Our nomenclature of "Hartree" and "Fock" diagrams is taken over from ref. [8]. There closed lir- and 2tt- 
exchange vacuum diagrams with two (resp. one) closed nucleon line have been called Hartree (resp. Fock) 
diagrams. The one-particle property S S pi n is obtained by opening one nucleon line of such diagrams. Conse- 
quently, we call here diagrams with one open and one closed nucleon line " Hartree" diagrams and diagrams with 
only one open nucleon line "Fock" diagrams. 



graph (right diagram in Fig. 1). In the static approximation the product of 7riV-interaction 
vertices a ■ {p — Pi) a • (p — Pi) — {p — Pi) 2 is spin- independent. A non-vanishing nuclear 
spin-orbit strength comes therefore only as a relativistic l/M 2 -correction. Isolating the ia x q 
factor from the product of fully relativistic pseudo-vector 7riV-interaction vertices, performing 
the 1/M-expansion, integrating over Fermi-spheres of radii k P:U = kf(l =F 5) 1 / 3 and expanding 
to linear order in the asymmetry parameter S, we get from the l7r-exchange Fock diagram: 



Here, we have introduced the dimensionless variables x = p/m n and u = kf/m n and the auxiliary 
function 

1 . l + (u + x) 2 

— In - - 

Ax 1 + (u — x) 



L(x, M ) = ^-ln ;;i" rJ ; 2 . (6) 



Next, we come to the evaluation of the four Hartree diagrams of iterated one-pion exchange 
shown in Fig. 2. We start with the left graph with one medium insertion, labeled (a). The 
relevant ia x q prefactor can be isolated already in the first step of the calculation from the 
product of 7riV-interaction vertices a ■ (I — q/2) a ■ (I + q/2) at the open nucleon line. For all 
remaining parts of that diagram one can then take the limit of homogeneous isospin-asymmetric 
nuclear matter (i.e. q — 0). After solving the inner <i 3 /-loop integral [6] and expanding the 
integrals over Fermi spheres of radii k p>n = kf(l =f5) 1 ^ 3 to linear order in 5, we find the following 
closed form expression for the isovector nuclear spin-orbit strength generated by the Hartree 
diagram (a): 

V^\p,k f ) = %~Stt{y arctan(u-x)-arctan(u+x)^+2x(^-u 2 +xA L(x,u)+2ux\ . (7) 

We continue with the calculation of the Hartree diagrams with two medium insertions, labeled 
(b), (c) and (d) in Fig. 2. In these diagrams, a (small) momentum transfer <f 12 with <f = g! + <f 2 
occurs at each medium insertion. The important prefactors ia x g} can again be isolated in the 
first step of the calculation from the product of 7riV-interaction vertices at the open nucleon line. 
After that the vectors cfj can be set to zero in all remaining components of these diagrams. When 
Fourier-transformed with density form factors to coordinate space each such momentum transfer 
qj, (j = 1,2) leads to the expression: (2ir)~ 6 f d 3 q J d 3 qj e^' r iqj &(qj) <&(q — qj) = f(r)Vf(r), 
with /(r) the density profile of weakly inhomogeneous nuclear matter. Consistent with the 
assumption of a weakly inhomogeneous nuclear matter (i.e. keeping only linear terms in small 
spatial gradients) we can make the replacement: f(r) V/(r) — > V/(r). The essential conclusion 
from this consideration is that for the calculation of S spin in weakly inhomogeneous nuclear 
matter each momentum transfer qj can be identified with q. With the help of this rule and 
certain techniques to reduce six-dimensional principal value integrals over the product of two 
Fermi-spheres, we end up after expanding to linear order in 8 with the following result for the 
Hartree diagram (b): 
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with the auxiliary functions s = xy + \fu z — x 2 + x' 2 y' 2 and s' — u ds/du. Throughout this work 
the momentum p is restricted to the interval < p < kf. Similarly, we find for the Hartree 



diagram (c): 



As' 



- xy Hsx[l + (xy + 2 ] 2 



s'(# + 2xy) 



u 



2< + (u 2 -£ 2 )ln 



l + (xy + 2 
and for the Hartree diagram (d): 

g\Mml f-,2e f 1 



u 



/\ 2i_ M + f 



(4s + s')w In 



(9) 



a; 2 — £ 2 y 2 



3(7 + 2<r 

1 + ff 2 



— 3 arctan er 



+ 



(y ln ^±44 _ 2,i [fa - ) + 5ff _ 6 arctan (j 



(l + a 2 )^ 



(10) 



with the auxiliary functions a = £y + Vw 2 — £ 2 + £ 2 ?/ 2 and a' = u da/du. The symbol -f\ dy in 
eq.(10) stands for a principal value integral. Note the similarity of the expressions eqs. (8,9,10) 
with the three-body isovector single-particle potentials written in eqs. (30, 31, 32) of ref.[7]. 
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Fig. 3: The iterated ln-exchange Fock graphs. These diagrams are labeled (e), (f), (g), (g) from 
left to right. 

Next, we come to the evaluation of the four iterated l7r-exchange Fock diagrams shown 
in Fig. 3. We start with diagram (e) with one medium insertion. The method to isolate the 
important proportionality factor q has been explained in section 3 of ref. [7]. Considering the 
isospin factor of diagram (e) and expanding the integrals over Fermi spheres of radii k p>n = &/(lT 
5) 1 / 3 to linear order in the asymmetry parameter 5, we end up with the following contribution 
to the isovector nuclear spin-orbit strength: 
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We continue with the computation of the Fock diagram (f) with a symmetrical arrangement 
of the two medium insertions. It leads to the following contribution to the isovector nuclear 
spin-orbit strength: 
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with the auxiliary functions G(x, u) = u(l + u 2 + a; 2 ) — [1 + (u + x) 2 } [1 + (u — x) 2 ]L(x, u) and 
t = xz + Vw 2 — a; 2 + x 2 2; 2 . Finally, we have to evaluate the last two topologically distinct Fock 
diagrams in Fig. 3. Since they contribute equally to the spin-dependent interaction energy E spin 
we have given both diagrams the same label (g). We find the following representation for their 
contribution to the isovector nuclear spin-orbit strength: 
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for the combined integrand of eqs. (21,23) in ref. [7] as well as the auxiliary function R = 
^J(l + x 2 — £ 2 ) 2 + 4£ 2 (1 — y 2 ). This completes the presentation of analytical results for the 
isovector nuclear spin-orbit strength Vi s (p, kf) generated by l7r-exchange and iterated l7r-exchange. 



3.1 Results 

For the numerical evaluation of the isovector nuclear spin-orbit strength Vi s (p, kf) we use con- 
sistently the same parameters as in our previous works [7, 8]. We choose the value g\ = 1-3 
for the nucleon axial vector coupling constant. The weak pion decay constant has the value 
f n = 92.4 MeV and M = 939 MeV and m n = 135 MeV are the masses of the nucleon and the 
(neutral) pion, respectively. 



Diagram 


lvr-Fock (a) (b) (c) (d) (e) (f) (g) 


V la (0,k fo ) 


-0.14 24.72 -1.49 -18.94 51.40 -42.35 -4.55 41.94 


Vi s (kfo, kf ) 


0.00 22.03 10.60 -10.64 25.83 -38.40 -2.68 27.04 



Table 1: Contributions of individual diagrams to the isovector nuclear spin-orbit strength 
Vis{p, kfo) atp = and at p = k f0 = 272.7 MeV. The units are MeVfm 2 . 

In the second row of Table 1, we present numerical values for the contributions of individual 
diagrams to the isovector nuclear spin-orbit strength Vi s (0, kfo) at the nuclear matter saturation 
density kfo = 272.7 MeV (as predicted by our calculation [8]). As expected the relativistic 
l/M 2 -correction from the l7r-exchange Fock graph is a very small effect. The contributions 
of individual iterated l7r-exchange diagrams are large and comparable in magnitude to the 
respective contributions to the isoscalar nuclear spin-orbit strength C/j s (0, kf ) (see Table 1 in 
ref.[7]). The basic reason for these large values is the large scale enhancement factor M (the 
nucleon mass) entering the iterated l7r-exchange. It stems from the energy denominator of 
such second-order diagrams which is a difference of small nucleon kinetic energies. Adding up 
the entries in the second row of Table 1 one gets VJ s (0, kfo) = 50.6MeVfm 2 , which is about 
1.4 times the isoscalar nuclear spin-orbit strength U(0,k ko ) = 35.1MeVfm 2 [7] generated by 
the same 27r-exchange diagrams. The predicted total sum is dominated by the contribution 
v is H \°, k fo) = 55.7MeVfm 2 of the four Hartree dia grams (a), (b), (c) and (d) (see Fig. 2). 
The same feature, namely the numerical suppression of the iterated l7r-exchange Fock diagrams 
against the Hartree diagrams, holds also for the isoscalar spin-orbit strength £4s(0, kfo) [7] as 
well as for the (spin-independent) nuclear mean- field studied in ref.[8]. The numerical entries 
(a)-(g) in Table 1 are all of the same order in the small momentum expansion. Individual 
cancelations give therefore no hint on the "convergence" of the chiral expansion of Vi s (p,kf). 
Note that the isovector nuclear spin-orbit strength VJ s (0, kfo) generated (almost completely) by 
iterated l7r-exchange is neither of relativistic nor of short range origin. It is in fact linearly 
proportional to the nucleon mass M and its inherent range is the pion Compton wavelength 
m^ 1 = 1.46 fm. The full pion Compton wavelength m" 1 and not half of it is suggested by 
comparing the expression V l ^ l \p,kf) in eq.(7) with the l7r-exchange single-particle potential 
U2(p, kf) in eq.(8) of ref.[8]. As a matter of fact the same non-polynomial functions L(x,u) and 
arctan(w ± x) occur in both cases. 

In Fig. 4, we show by the full line the dependence of the calculated isovector nuclear spin- 
orbit strength Vi s (0,kf) on the nucleon density p = 2k 3 f/3n 2 . One observes in the region 
p < 0.4 fm -3 an approximate linear growth of VJ s (0, kf) with the density p as it is typical 
for relativistic mean-field models. In a finite nucleus the spin-orbit force acts mainly on the 
surface where the density gradients are largest and the density has dropped to about half of 
the central density. The replacement f(r) V/(r) — > |V/(r) (instead of f(r) V/(r) — > V/(r) 
valid for weakly inhomogeneous nuclear matter) describes then more realistically the situation 
for a finite nucleus. The dashed line in Fig. 4 shows the isovector nuclear spin-orbit strength 
which results if the contributions from the diagrams with two medium insertion (b), (c), (d), (f) 
and (g) are weighted with a factor 1/2. This different weighting reduces of the total isovector 
spin-orbit strength by about a factor 3. 
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Fig. 4-' The full line shows the isovector nuclear spin-orbit strength Vi s (0,kf) at zero nucleon 
momentum (p = 0) versus the nucleon density p = 2A;|/37r 2 . If the diagrams with two medium 
insertions are weighted with a factor 1/2 the dashed line results. 
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Fig. 5: The momentum dependence of the isovector nuclear spin-orbit strength Vi s (p,kfo) at 
saturation density kfo = 272.7 Me V. The full (dashed) line corresponds to a weighting of the 
diagrams with two medium insertions with a factor 1 (1/2). 

In Fig. 5, we show by the full line the dependence of the calculated isovector nuclear spin- 
orbit strength Vi s (p, kfo) at saturation density kjo = 272.7 MeV on the nucleon momentum p for 
< p < kfo- One observes a moderate p-dependence. In contrast to the isoscalar nuclear spin- 
orbit strength Ui s (p, kf ) shown in Fig. 7 of ref.[7] there is no sign change in the present case. 



The dashed line in Fig. 5 corresponds to the weighting of diagrams with two medium insertions 
with a factor 1/2. It shows the same qualitative p-dependence as the full line in Fig. 7. The 
entries in the third row of Table 1 indicate how the contributions from individual diagrams vary 
with the nucleon momentum from p = to p = kfo. 

The terms from iterated l7r-exchange written down in eqs.(7-14) are the unique leading order 
contributions in the small momentum expansion. They are of second order in small momenta 
as signaled by the common prefactor g\Mm1/(2nf n ) 4: . Nuclear structure calculations which 
use the here calculated isovector nuclear spin-orbit strength Vi s (p,kf) as input would be most 
useful in order to clarify the role of the spin-orbit interaction generated by 27r-exchange in a 
finite nucleus. 

4 Relativistic corrections to the isoscalar nuclear spin- 
orbit strength 

In this section we calculate several higher order corrections to the isoscalar nuclear spin-orbit 
strength Ui s (p, kf) generated by 27r-exchange in order to learn about the convergence of the un- 
derlying small momentum expansion. All these contributions scale like 1/M with nucleon mass 
M and can therefore be regarded as truly relativistic spin-orbit effects. Since our calculation 
is a strictly perturbative one 7riV-interaction vertices in-medium are the same as in vacuum. 
The relativistic spin-orbit effects depend therefore reciprocally on the free nucleon mass M and 
not on a self-consistently determined effective nucleon mass M* as it is the case in relativistic 
mean-field models [2]. 




Fig. 6: Some diagrams representing the isoscalar nuclear spin-orbit interaction arising from 
irreducible two-pion exchange. The first two (last two) diagrams are constructed from the 2n- 
exchange triangle (crossed box) graph of elastic NN-scattering. 

In Fig. 6 we show some diagrams which generate the isoscalar nuclear spin-orbit strength 
related to irreducible two-pion exchange. By opening them at the medium insertion (symbolized 
by the double line) ones see that the first two Fock diagrams stem from the 27r-exchange triangle 
graph of elastic NN-scattering, whereas the last two Hartree and Fock diagrams are constructed 
from the 27r-exchange crossed box graph of elastic NN-scattering. We have explicitly calculated 
the isoscalar nuclear spin-orbit strength Ui s (p,kf) generated by these diagrams using the tech- 
niques outlined in section 2. In an intermediate step the resulting isoscalar nuclear spin-orbit 
strength Ui s (p, kf) can be written as an integral of the spin-orbit NN-amplitudes of the triangle 
and crossed box graph [6] over a Fermi sphere of radius kf. This expression can be further 
reduced to a one-dimensional integral of the form: 

Ut\p,kf) = j£'%#J-(p+g) 2 ^ ( 15 ) 



with Vso(q) and Wso{o) the isoscalar and isovector spin-orbit NN- amplitudes. The term 
2^50(0) in eq.(15) belongs to Hartree type diagrams (with one closed and one open nucleon 
line) while the term Vso(<?) + 3Wso(<l) summarizes the contributions from Fock type diagram 
(having just one open nucleon line). The spin-orbit NN-amplitudes related to irreducible 2n- 
exchange have been calculated in section 4.2 of ref. [6] using dimensional regularization. If we 
use instead cut-off regularization the linear divergences of the irreducible 27r-exchange diagrams 
become visible. We find: 
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which agrees up to the additive constants linear in the momentum cut-off A with the expressions 
in eqs. (22,23) of ref. [6]. Note that eqs.(16,17) include also the irreducible contributions of the 
27r-exchange planar box graph (not shown in Fig. 6). Inserting eqs. (16, 17) into the "master 
formula" eq.(15) one gets for isoscalar nuclear spin-orbit strength arising from irreducible two- 
pion exchange: 
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At zero nucleon momentum (p = 0) this simplifies to: 
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Note the reciprocal dependence on the nucleon mass M which indicates that this nuclear spin- 
orbit strength is a relativistic effect. 
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Fig. 7: The full line shows the isoscalar nuclear spin-orbit strength Uig n \p, kf ) related to irre- 
ducible "In-exchange versus the nucleon momentum p . The dashed line gives the combined result 
of the three-body diagrams in Fig. 8 and the dashed-dotted line corresponds to the relativistic 
1/M -correction to the iterated ln-exchange Hartree diagram (a) in Fig. 2. 

As in our previous work [8] we choose the value A ~ 0.65 GeV for the momentum cutoff 
A. This gives at saturation density kf = 272.7 MeV an isoscalar nuclear spin-orbit strength 

due to irreducible 27r-exchange of U^ n) (0, k f0 ) = (16.20 - 9.27) MeVfm 2 = 6.93MeVfm 2 . This 
combination of A-dependent an finite terms is less than 20% of the total isoscalar spin-orbit 
strength coming from iterated l7r-exchange. Such a reduction factor is in accordance with the 
expectations from chiral power counting since the first/second term of 11^(0, kf) in eq.(19) 
is suppressed by one/two powers of small momenta in comparison to the results from iterated 
l7r-exchange. The full line in Fig. 7 shows the momentum dependence of the isoscalar nuclear 
spin-orbit strength U^\p, kf ). One observes a very weak fall-off with momentum from p = 
to p = k f0 = 272.7 MeV. 




Fig. 8: Three-body force Fock diagrams constructed from the (isovector) Weinberg- Tomozawa 
ttitNN -contact vertex. Analogous Hartree diagrams which have an isospin factor zero in isospin- 
symmetric nuclear matter are not shown. 



By putting a second medium insertion to the intermediate nucleon propagator in the first 
two diagrams in Fig. 6 the first two graphs shown in Fig. 8 result. By opening them at the two 
medium insertions one sees that these graphs can be classified as three-body diagrams. With 
employment of one (isovector) Weinberg- Tomozawa 7T7riViV-contact vertex one can construct a 
tree-level 27r-exchange three-body diagram (NNN— >NNN). The three Fock diagrams in Fig. 8 
result then from this three-body diagram by closing two nucleon lines in all possible ways. 2 
Let us begin with the computation of the right Fock diagram in Fig. 8 having a symmetrical 
arrangement of the two medium insertions. The 7T7riViV-contact vertex (symbolized by the heavy 
dot) contributes for the relevant kinematics first in form of a relativistic 1/M-correction. After 
isolating the proportionality factor iff x q from all spin-dependent terms one is left with an 
integral over the product of two Fermi spheres of radius kf. The latter can be solved in closed 
form and one gets the following contribution from the third diagram in Fig. 8 to the isoscalar 
nuclear spin-orbit strength: 
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with the auxiliary function G(x,u) defined after eq.(12). The computation of the first two 
(topologically distinct) diagrams in Fig. 8 is done in a similar fashion. In this case one non- 
elementary integration is left over from the six-dimensional integral over the product of two 
Fermi spheres of radius kf. We find the following expression for the isoscalar nuclear spin-orbit 
strength: 
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to which the first and second diagram in Fig. 8 have contributed equally. Numerical evaluation 
of eqs. (20,21) gives at nuclear matter saturation density kf = 272.7 MeV and zero nucleon 

2 Wc discard the Hartrcc-type diagrams resulting from this closing procedure since they have an isospin factor 
zero in isospin-symmetric nuclear matter. 



momentum U^ ] {0, k f0 ) = (-0.20 - 1.10) MeVfm 2 = -1.30MeVfm 2 . This is a very small effect 
compared to iterated l7r-exchange. The momentum dependence of the combined result from 
eqs. (20,21) for uj®(p, kf ) is shown by the dashed line in Fig. 7. One observes a weak increase 
with momentum from p = to p = kfo- In the heavy baryon framework the relativistic 1/M- 
corrections calculated here are generated by (nucleon-momentum dependent) vertices from the 
second order chiral Lagrangian C^ N (proportional to 1/M). The other 7T7riViV-contact vertices 
proportional to the low-energy constants c 12 ,3,4 give rise to spin-orbit terms first one order higher 
in the small momentum expansion. An order of magnitude estimate of their effect on U[ s (p, kf) 
can be obtained by inserting the spin-orbit NN-amplitudes Wso(o) and Vso(o) in eqs. (13,15,16) 
of ref. [13] (which do not include regularization dependent polynomial pieces) into the "master 
formula" eq.(15). 

There are of course many more contributions to the isoscalar nuclear spin-orbit strength at 
fourth order in small momenta (having the common prefactor m 4 /M). In order to get an idea 
about the typical size of other such effects we consider now the first relativistic 1/M-corrections 
to the dominant iterated l7r-exchange Hartree diagrams in Fig.2. In order to obtain conveniently 
all 1/M-corrections from these diagrams one starts with fully relativistic interaction vertices and 
fully relativistic pion- and nucleon propagators sandwiched between Dirac spinors. Next one 
performs an expansion in powers of 1/M to the appropriate order and solves the remaining 
(loop)-integrals. In the case of diagram (a) in Fig. 2 with one medium insertion one ends up 
with the following 1/M-correction to the isoscalar nuclear spin-orbit strength: 
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In this expression we have subtracted the irreducible part g\k 3 f (6A — 57rm 7r )/(1287r 4 M/^) which 
is already accounted for in eqs. (15,18) by the spin-orbit NN-amplitude Vso(0) derived from 
the planar 27r-exchange box graph. Inserting the cut-off scale of A ~ 0.65 GeV one gets 
f// s a) (0, k f0 )\ M -i = (-0.85 + 3.41) MeVfm 2 = 2.56MeVfm 2 which is a small -5% correction 
to the leading order contribution [7^(0, &/o)|m = —52.96 MeVfm 2 [7]. The momentum depen- 
dence of Uig\p, kf )\ M -i is shown by the dashed-dotted line in Fig. 7. One observes a moderate 
increase with momentum from p = to p = A; jo- 
in order to avoid very lengthy analytical expressions, we present here only the first rela- 
tivistic 1/M-corrections to the Hartree diagrams (b) and (c) in Fig. 2 evaluated at zero nucleon 
momentum (p = 0). In the case of diagram (b) one finds: 
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The resulting numerical value 11^(0, fc/o)^- 1 — —1.06 MeVfm 2 is again a small —3% correction 
to the leading term U^(0, k f0 )\ M = 37.35 MeVfm 2 [7]. The analo srous result for the Hartree 



diagram (c) reads: 
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which leads to an even smaller numerical value of 11^(0, /c/o^m- 1 — 0.10 MeVfm 2 . 

In summary, we have calculated here several relativistic 1/M-corrections to the isoscalar 
nuclear spin-orbit strength U[ s (p, kf) generated by chiral 27r-exchange. As expected from chiral 
power counting these higher order terms are significantly smaller than the unique leading contri- 
bution from iterated l7r-exchange [7]. Its role in finite nuclei should be clarified by implementing 
it into nuclear structure calculations. 
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